Swift J164449. 3+573451 event: generation in the collapsing star cluster? 
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We discuss the multiband energy release in a model of a collapsing galactic nucleus, and we try to 
interpret the unique super-long cosmic gamma-ray event Swift J164449. 3+573451 (GRB 110328A by 
early classification) in this scenario. Neutron stars and stellar-mass black holes can form evolutionary 
a compact self-gravitating subsystem in the galactic center. Collisions and merges of these stellar 
remnants during an avalanche contraction and collapse of the cluster core can produce powerful 
events in different bands due to several mechanisms. Collisions of neutron stars and stellar-mass 
black holes can generate gamma-ray bursts (GRBs) similar to the ordinary models of short GRB 
origin. The bright peaks during the first two days may also be a consequence of multiple matter 
supply (due to matter release in the collisions) and accretion onto the forming supermassive black 
hole. Numerous smaller peaks and later quasi-steady radiation can arise from gravitational lensing, 
late accretion of gas onto the supermassive black hole, and from particle acceleration by shock waves. 
Even if this model will not reproduce exactly all the Swift J164449. 3+573451 properties in future 
observations, such collapses of galactic nuclei can be available for detection in other events. 

PACS numbers: 98.54.Cm, 98.70.Rz 
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I. INTRODUCTION 

On March 28, 2011, the Swift's Burst Alert Telescope 
detected the unusual super-long gamma-ray event Swift 
J164449. 3+573451 from the center of a star-forming 
galaxy at the redshift z = 0.3534+0.0002 Q. The burst 
sustained huge activity for a few days, and it keeps a 
softer but steady activity for a few months. On the first 
days, the event had an extremely variable light-curve 
with several high peaks, numerous smaller spikes, and 
a quasi-continuum. The peak luminosity of this burst 
reached a level of ~ 5 x 10 erg s~ x . It is doubtful that 
this event can be attributed to the class of superlong 
(with a duration > 500 s) GRBs 0, instead it seems to 
be a completely different phenomenon. 

As a possible interpretation the of Swift 
J164449. 3+573451 the tidal disruption flare model 
of 3 was adopted in 0-Q|- I* 1 this model a star 
was destructed tidally by the supermassive black hole 
(SMBH) with a mass of ~ 1O 7 M in the galactic centre, 
the debris was accreted onto the SMBH, and it supposed 
that a jet was produced in a direction almost precisely to 
the Earth. The later configuration is similar to a small 
blazar |8j. The accidental jet direction to the Earth in 
Swift J164449. 3+573451 for producing the huge variable 
activity. 

Here we propose an alternative explanation, based on 
the model of gravitational collapse of a central compact 
cluster and multiple mergers of neutron stars (NSs) and 
black holes (BHs) there. Dense star clusters are a typical 
structural element of galactic nuclei. Compact subclus- 
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ters of stellar remnants can form in the course of evolu- 
tion due to mass segregationsinking of the most massive 
stars to the center and their explosions as supernovae 0- 
ITTj ]. We hereafter refer to this inner subsystem of com- 
pact objects as a "cluster". The cluster evolves dynam- 
ically due to two-body relaxation: the central density 
and velocity dispersion grow, fast objects "evaporate" 
from the central core, forming a relaxed isothermal halo 
around the central core. The central part (core) of this 
cluster approaches a relativistic stage, and the possibil- 
ity of gravitational collapse is limited by relaxation time 

There are several possible mechanisms of gamma ray 
generation and energy release in other bands during and 
after the core collapse: 

1. Collisions of neutron stars and stellar- mass black 
holes can directly generate gamma-ray bursts, as is the 
case in the ordinary model of short GRBs. Before the 
collision, two remnants are captured into a close orbit, 
and jet formation in this system can proceed similar to 
the numerical calculations of [133- The fastest variability 
time-scale of Swift J164449.3+573451 is - 100 s. So the 
separate GRBs could explain the high peaks in the light- 
curve of Swift J164449. 3+573451 only if there is some 
mechanism of the peaks broadening as compared with 
usual short GRBs having durations < 2 s. This can be 
due to a high gas density in the collapsing galactic nu- 
cleus. The prolonged emission lasts ~ 10 4 s in the model 
of ordinary and compacts stars collisions [l5[ , [l6| . In the 
case of NS-NS and NS-BH collisions the extended emis- 
sion can be sustained for ~ 100 s, as it was shown in the 
two jet model of [13]. Possibly, direct gamma pulses were 
not observed because the Earth was not located in the 
collimated jet directions, in this case the followingmech- 
anisms can be responsible for peak generation. 

2. Each collision of compact objects (NS-BH) provides 
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matter release in the form of debris, fireballs and jets. 
This is an alternative source of gas accretion in addition 
to star tidal disruption 0-0]; but the matter of a neu- 
tron star produces denser and more compact (initially) 
ejection. Collisional destructions of several neutron stars 
provide a larger power then a single ordinary star. Just 
after the SMBH formation, several (from different col- 
lisions) mutually fighting accretion disks or one steady 
disk (along the angular momentum of SMBH) can form. 
The bright peaks in the Swift J164449.3+573451 light 
curve during the first two days may therefore be a conse- 
quence of multiple episodes of violent gas accretion onto 
the forming SMBH. Each NS collision and destruction 
is responsible for a separate peak in the light curve of 
Swift J164449. 3+573451. Then the accretion disk forms, 
and different radiative processes develop in a similar way 
to the ordinary scenario with tidal destruction of a star. 
Therefore many properties of the tidal destruction model 
could be relevant to our model with only minor modifi- 
cations. 

3. Numerous smaller peaks can arise due to gravita- 
tional lensing of the signals (from mergers and accretion) 
of the forming SMBH and compact objects in the cluster. 
Later lensings are similar to "mesolensings" [l8| and can 
produce numerous subpeaks in the light curve. 

4. Smooth flashes and quasisteady radiation can arise 
from gas accretion onto the SMBH at the time of its 
formation and rapid growth (for a few first days) and 
later. 

5. The galactic nucleus as a whole will be source of 
IR radiation because of multiple fireball formation in- 
side the gas-dust envelope. The debris of disrupted nor- 
mal stars form the envelope in which the NS/BH clus- 
ter is submerged. Really, the X-ray spectrum of Swift 
J164449. 3+573451 shows strong absorption which can 
witness a wide envelope around the central engine. Nev- 
ertheless, to explain Swift J164449. 3+573451 one must 
suppose that this gas envelope is optically thin for gamma 
radiation. Acceleration of particles by multiple shock 
fronts of the fireballs is also possible. 

6. This model also predicts huge multiple bursts of 
gravitational waves . Similar events will be detectable 
in the near future by new gravitational wave telescopes. 

The first of the effects listed above is similar to the 
GRB recurrence scenario which has been elaborated in 
[20I - I22I ] - This model associates the recurrent and mul- 
tiple GRBs with the final dissipative stage of central 
star clusters in galactic nuclei, which precedes SMBH 
formation. Accidental gravitational captures produce 
short-lived pairs of compact objects. According to the 
most popular scenarios, short GRBs (with a duration of 
< 2 s)are generated during coalescence of two NSs or a 
NS and a BH in binary systems [23|, [24{. It has been 
shown in recent calculations |l3j that mergers really pro- 
duce short GRBs. Similar but multiple GRBs can be 
generated during the clusters core gravitational collapse 
just before the SMBH birth. This model could also ex- 
plain the unusual double burst GRB 110709B [3. 



II. EVOLUTION OF COMPACT STAR 
CLUSTER 



According to observations, about 50-80% of galaxies 
with low and intermediate luminosities contain the com- 
pact central star clusters [25j, as well as our Milky Way 



26] . Most massive stars of these star clusters inevitably 
sank down to the dynamical centre and exploded as su- 
pernova long ago at the initial epoch of the star clus- 
ter evolution. Therefore, the formation of the compact 
subsystems of NSs and stellar mass BHs must be a com- 
mon process in the Universe [§], [l(|, [HI- Let us con- 
sider the cluster of NSs and BHs, formed at the time 
ti <C to, 35, where £0.35 — 10 Gyr corresponds to the red 
shift z = 0.35. The subscript "i" marks the quantities 
at the moment ti. During the dynamical evolution of 
the cluster's core its mass M c , number of compact ob- 
jects N c = M c /m, velocity dispersion v, virial radius 
R c = GN c m/2v 2 and the masses of separate objects m 
are all evolving. We use the homological approximation 
for the mean properties of the core. In this approxima- 
tion the behaviour of the core obeys the system of equa- 
tion, obtained in [l9[ . The global parameters of the core 
change due to objects dynamical evaporation from the 
core and due to energy loss for gravitational radiation, 
while the rate of the above processes and the inner struc- 
ture of the cluster are determined by the two-body relax- 
ation and mergers of the compact objects. We rewrite the 
system of equations from [191 ] through the new variables 
x = v 2 /c 2 in the following form: 



x/x — (ct2 — a\)t r 
M/M = -a 2 t~ 1 
rh/m = t~ 



+ 7«-y5 



-1 

"cap' 



(1) 



where i cap takes into account the dissipative processes 
of mergers of compact objects and gravitational waves 
emission, t~JL = a ca/P v/y/2, and t r is the two-body re- 
laxation time (characteristic time-scale for stellar cluster 
dynamical evolution) 
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where A = ln(0.4jV) is the Coulomb logarithm and we 
put A ~ 15 hereafter, i. e. A ~ const with logarith- 
mic accuracy, n is the compact objects number density. 
According to [27J, the capture cross-section for the two 
objects with masses m\ and m,2 is 
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= m in this Section. 



(3) 



The coef- 

-3 



and we put m\ 

ficients a x = 8.72 x 10~ 4 and a 2 = 1.24 x 10"" were 
obtained in the Fokker-Plank model [19|. This formal- 
ism describes the contraction of the clusters' core, while 
the evaporated objects form the isothermal halo around 
the core with density profile p oc r~ 2 . The halo can 
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be involved into the avalanche contraction after the core 
reaches a relativistic stage, as will be discussed in the 
next section. 

Most of the cluster lifetime is non-dissipative (one can 
neglect the t cap terms). From the solution of (pj in this 
case we have: 

M c = M i [l-(t-t i )/(Kt ri )] l '\R c = R4l-(t-U)/{Kt ri )p, 

(4) 

where v\ — 2a 2 /(7c(2 — 3ai), v 2 = 2(2a 2 — cn)/(7a2 — 
3ai), K = 2/(7a 2 — 3ai). The quantity t e = Kt r .i ~ 
330£ r .i approximately equals to the duration of the clus- 
ter evolution up to the time of the full evaporation or 
gravitation collapse into a SMBH, because the final dis- 
sipative stage is rather short. 

The solution of the exact Equations (TTJ has the form 



m (x) = m t (x 5 / 7 ' + x%')/(x°" + <'') 



5/7s 



„ 5 /7 , 5/7, 



M c (x) = Mi 



*wj+* b/7 ) 

x^{x*£ + xf) 



St»2->l) 



where 



10A(a 2 -ai)/(7v^3) 



7/5 



5.3 x 10" 



(5) 
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The cluster of NS had formed long before the dissipative 
stage. Near the dissipative stage x a;,. Putting m = 
2m i in ([5]), we obtain x — Xdis- Therefore, the dissipative 
stage begins at v/c ~ 0.023, i. e. v ~ 6.9 x 10 3 km s _1 , 
independently from the initial parameters of the cluster. 
Duration of the dissipative stage 



of BH have time to grow up to ~ rrii(xf /xdi s ) 5 ^ 7 ~ 42wi, 
according to ((SJ), although the validity of ([SJ is doubtful 
near the moment of the collapse. The homological model 
gives a useful qualitative description, but numerical re- 
sults must be accepted with some caution. This picture 
of the dynamical evolution is somehow complicated in the 
presence of an angular momentum, but it is not impor- 
tant in the central relaxed part of the cluster. A serious 
modification is required for a cluster with a pre-existing 
massive BH at the center. In this case the central BH 
grows slowly instead of a fast avalanche collapse of the 
cluster. 

In plenty of works, the dynamical evolution of globular 
clusters was considered in the "evaporation model" (see, 
e - g-> [HI)- This model supposes that an object leaves the 
core with zero kinetic energy, so that the total energy of 
the evolving core is conserved E c oc M^/R c « const. 
The rate of evaporation was calculated from the high- 
velocity tail of Maxwellian distribution. Let us compare 
the evaporation model with the considered above model 
of In the evaporation model a\ = and a 2 = 

7.4 x 10 -3 , and the duration of core evolution till the 
collapse t e ~ 2/(7a 2 ) ~ 38.6t r ,i is several times less in 
comparison with the above calculations. It means that 
the core more easily reaches the relativistic stage and its 
initial mass Mi could be smaller. For this model Xdis — 

3.5 x 10~ 2 (v ~ 0.19c), r dis ~ 3.6 yr, Xi ~ 2 x 10~ 7 
(vi ~ 140 km s" 1 ) and M min ~ 3 x 10 6 M Q . 

One can roughly suppose that the epoch of collapses 
lasts to ~ 10 Gyr, then the total rate of such collapses in 
the observable Universe is estimated to be 



N h - —(ctofn otn 1 w 0.1 { 

3 V 01 9 V 10- 2 Mp C - 3 



( — ff^yr" 1 , (ii) 
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Tdi: 



(Xdis) 



( M c 



\10 a MqJ \2.8M, 



yr- (8) 



Collapse occurs at x ~ xj ~ 0.1 1!) . If we require the 
collapse at z = 0.35, i.e. Kt r ^ — to. 35 — 9.8 x 10 9 yr, we 
get 



where n g is the number density of structured galaxies 
with central clusters. It seems natural that the collapse 
predicted by this model could be seen for several years 
of observations. Moreover, the clusters core just before 
the collapse can be a promising source of multiple gravi- 
tational wave bursts 1191. 



9 x 10 



-7 ( Mi 



V 10 7 M P 



4/3 



1.4M Q 



-2/3 



(9) 



and it corresponds to Vi ~ 280 km s _1 . If we require 
additionally that at the beginning of the collapse the core 
hasn't evaporated, M c (xf) > 2mj, then we have from ([6]) 
and © 



Mi > M n 



3 x 10 7 M G 



(10) 



or Ni > 2 x 10 7 if mi — 1.4M . The boundary mass 
■^min sets the minimum for the initial mass of the clus- 
ter in the scenario under consideration. Mass segregation 
accelerates the evolution as compared with this homolog- 
ical model. Therefore, the real M m ; n may be somehow 
smaller. We put in the following estimations M m j n ~ 
10 7 M Q . Note that till the moment i*~ 0.1 the masses m 



III. AVALANCHE CONTRACTION AND 
GRAVITATIONAL COLLAPSE OF THE 
CENTRAL CORE 

Now we consider the final highly dissipative stage of 
the cluster evolution. We suppose that the central clus- 
ter consists of stellar-mass BHs with an admixture of 
non-coalesced NSs. The BH masses in the cluster grow 
through mergers fl9j | . Suppose that their masses imme- 
diately before the core's collapse are torh ~ lOm^, where 
nii ~ IAMq. In the central core ttibh ~ 42mi according 
to the rough estimation of the previous section, but we 
put mBH ~ 10mj as mean over the cluster. We put also 
that the cluster has the total mass M (mainly in the form 
of BHs) and the NSs constitute a fraction /ns -C 1 of M. 
The dynamical evolution of clusters is accompanied by 
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the secular growth of the central velocity dispersion or, 
equivalently, by the growth of the central gravitational 
potential. Just near the collapse stage the cluster has a 
dense core with mass M c , radius R c ~ 3i? ff , c (last stable 
orbit), where R g _ c — 2GM c /c 2 , and velocity dispersion 
~ 0.3c. All objects with a small angular momentum will 
fall onto the core without returning to the outer part 
of the cluster. The avalanche occurs due to the pres- 
ence of quasielliptical orbits that connect different layers 
of compact objects 28]. Objects from the loss-cone are 
captured by the central core, and the central region col- 
lapses into the SMBH. Numerical calculations [29j] shows 
that the several percent of the total cluster mass will 
collapse onto the forming SMBH during the several dy- 
namical times idyn = Rd/v ~ duration of the avalanche 
contraction, where R c \ is the radius of the whole cluster 
(halo around the core). 

In view of /ns <C 1, the NS-BH coalescences are 
more probable than the NS-NS ones. At the mass ra- 
tio ttjbh ~ lOmj, the NS does not fall into the BH as 
a whole but is disrupted by tidal forces and produces a 
fireball and GRB. The rate of NS-BH coalescences in the 
core is N c = cr cap uJVnsIbHj where the gravitational cap- 
ture cross section for objects with different mass is given 
by (O, A^ns = fNsMc/m is the number of NSs in the 
core, and tt-bh is the BH number density in the core. For 
the particular parameters M — 10 M Q and v = 0.05c we 
have R c \ ~ 10~ 4 pc and idyn — 2.3 days. The number of 
NS-BH mergers in the core during this time is 



host galaxy of Swift J164449. 3+573451 didn't show activ- 
ity long before the event. The first activity, which could 
be associated with the beginning of the collapse and the 
fast collisions, was found in the Burst Alert Telescope's 
data at 3 days before the main event, and the flux of the 
precursor was ~ 7% of the first main peak Q . 

Just as the central relativistic region had appeared it 
began to capture objects which flew through this region 
inside the "loss-cone". For the object at the radial dis- 
tance r from the core's centre the loss-cone is character- 
ized by the angle ~ b/r, where b 2 ~ R g ^ c r p (c/v) 2 and 
the value of r p is the same as in the process of gravita- 
tional capture of two objects (27| with masses mi — tobh 
and 7Ti2 = M c m^: 

85ttG 7 / 2 tobhM c 5/2 

12V2c 5 v 2 ' [ ' 

For the calculations we use the isothermal density profile 
of the relaxed cluster p(r) — M / (AirR c \r 2 ) till the core 
radius R c and the isotropic velocity distribution. The 
number of BHs inside the "loss-cone" is estimated to be 



N c t dvn ~ 7 x 



M 



2 x 10" 5 J V 10 7 M ra / V 5 x 10 5 M 



M r 



(12) 

The rate of NS-BH collisions in the entire cluster before 
the core collapse: 



N ~ 2.5xl0~ 



/ns 



2 x 10- 5 



M 



10 7 M 



( 0.05c) 



31/7 



yr 



(13) 

i.e., the cluster doesn't show the fast recurrence before 
the collapse. Multiple short GRBs can be generated in 
the process of avalanche- like contraction at the cluster 
center [Hj], HI- Therefore, we have the several NS-BH 
collisions and, accordingly, several GRBs-likc events over 
the time of several days. Similar temporal characteris- 
tics were required in [221 ] to explain the multiple GRB 
recorded on October 27, 1996. 

In approach to the collapse moment t co \\ the rate of 
compact objects collisions (expression like ([131) ) rises 
abruptly, and the asymptotic solution of (JlJ near t co ii 
gives N c (x [1 — (t — t co \\)/T*]~ 2 , where the characteristic 
time-scale is 



0.007 x 



7->2 



CTgX 



31/14 



0.8 days, 



(14) 



where r g = 2Gtobh/c 2 . Therefore the stage of frequent 
NS-BH collisions lasts ~ one day. It explains why the 
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(16) 

For the fiducial parameters, used above, we have numer- 
ically N cone ~ 3 x 10 5 . Therefore, the mass inside the 
loss-cone is sufficient for our model. 

After the depletion of orbits inside the loss-cone dur- 
ing the several dynamical times the objects fall onto 
the SMBH only due to the slow diffusion of their or- 
bits through the loss-cone boundary with the rate Adif — 
N/(Xt r ), where t r is given by and A ~ ln(i? cl /3i? g . c ). 
For the parameters in use the loss-cone is empty, be- 
cause Adifidyn < A con , and at the post-collapse stage the 
SMBH's mass grows only slowly with the rate AjifmBH- 



IV. RADIATION FROM THE COLLAPSING 
NUCLEUS 

Let us consider different ways of energy release during 
and after the avalanche contraction and SMBH formation 
in a cluster of compact objects. 



A. Direct gamma-radiation 

NS-BH collisions can produce the short GRBs. If the 
gamma-radiation is collimated into the f2 7 solid angle, 
the number of the observed GRBs will be fi 7 /4-7r fraction 
of the number of NS-BH collisions (fl~2j). Let us denote 
the energy release at each merge by E — 10 52 E 52 erg. 
Then the mean power of the source during the active 
stage tdyn in the case of isotropic emission is 



L 



GRB 



= E N, 



2.3 x 10 47 £ 5 2 



/ns 

io- 4 



erg s 



(17) 
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This value coincides with the average luminosity of Swift 
J164449. 3+573451. For the larger / NS ~ 10" 3 - 10~ 2 
the N c is also greater, and one can achieve even the 
observed peak luminosity ~ 5 x 10 48 erg s _1 of Swift 
J164449. 3+573451. Alternatively, the peak signal can 
come from the single NS-NS or NS-BH merge events with 
beamed radiation [l3j]. In this case the high luminosity 
is the result of collimation. But one must require some 
mechanism of the pulses broadening from the typical 
~ 2 s duration of the short GRBs till the minimum time- 
variability scale - 100 s of the Swift J164449.3+573451. 
The mechanisms for the prolonged or extended emission 
after the first short bursts were discussed in 15 1, [l6j], 
[13. But the light-curve of the Swift J164449.3+573451 
has no short ~ 1 s pulses. A simple explanation of their 
absence is that we are not at one of the jets directions, 
and therefore another mechanism must be responsible for 
radiation with ~ 100 s variability. 



B. Radiation from the debris accretion 

Collisions between NSs and BHs provide a release of 
matter which can be accreted as in the star tidal dis- 
ruption models (343- Ejections from NSs could be more 
dense and compact, so they could provide a larger power 
as compared with a single ordinary star. Tidal disrup- 
tion of a non-ordinary star (a white dwarf) tidal disrup- 
tion by the SMBH with mass < 1O 5 M0 as a model of 
Swift J1644+57 was proposed in [3(| and the short time- 
scales in the light curve were explained. Here we pro- 
pose an even more extreme model with the NSs multiple 
destructions. The difference is that NSs destructed not 
by the tidal forces of central SMBH, but in the collisions 
with stellar mass BHs in the cluster. In this case the char- 
acteristic time-scales ~ 100 s are related to the SMBH 
horizon light crossing time 2GM/cr as in the model with 
the star tidal destruction. 

Note that the orbits of NSs in the cluster are dis- 
tributed almost isotropically. This means that the de- 
bris from the collisions fall onto the central SMBH from 
different directions. If each collision produces a separate 
short-lived accretion disk around the SMBH, the corre- 
sponding jets from different disks will also be pointed 
isotropically. We can reproduce the configuration of a 
small blazar with an accidental jet direction to the Earth 
in this case only by requiring that the actual number 
of collisions was larger then the number of the observed 
peaks by a factor of 47r/f2j, where Oj is the solid angle 
of the jet collimation. 

Another possibility states the existence of a steady ac- 
cretion disk around the SMBH that supports a single jet 
in the blazar configuration with a matter feed from the 
NS-BH collisions. In this model the jet direction is not 
related to the orientation of the debris orbits, but can 
be determined, for example, by the SMBHspin that had 
originated from the common rotation of the parent clus- 
ter. Radio and other signals can be generated as in the 



model with star tidal destruction, and we will not repeat 
the corresponding calculations here. 



C. Gravitational lensings 

Let us first consider the lensing of the point source 
on the point lens in the case the both are inside the 
same cluster at the distances from the Earth D s and 
Di respectively, so that D s « Di rj 1 Gpc (galaxy at 
z = 0.35), and D s — D\ is of the order of the cluster 
radius R c \ ~ 10 -4 pc. This configuration is similar to 
the "mesolensings" in globular clusters fl8| . When the 
Einstein radius of the lens is 



R E = (AGMi/c 2 ) 1 ' 2 [D,(l - Dt/D.)] 



1/2 



4 x 10" 



V lOMp 



1/2 



3 x 10 14 cm 



1/2 



(18) 



where Mi is the lens' mass. If one uses the r~ 2 den- 
sity profile of the cluster and suppose that the collision 
occurred in the core, then the integration gives the num- 
ber of lensing BHs inside the Einstein radius at the level 
Ni ~ 5 x 10 -4 . A similar small chance is in lensing on 
the central supermassive BH because the large mass is 
compensated by a small distance from the source to the 
lens (collisions are most probable in the core, near the 
central BH). This means that lensing of a point source 
or jet on BHs is unlikely. 

The large fireball lenses easily, because the radius of 
the cluster R c \ ~ 10~ 4 pc is of the order of the size i? sw 
of shock waves in GRBs where some part of the gamma- 
radiation are generated. First of all, the fireball scat- 
tered on the forming central SMBH. The hot medium 
flew around the SMBH and was swallowed inside it. Let 
the observer and an NS-NS or NS-BH merging be on op- 
posite sides of the SMBH. The surface of the fireball is not 
spherical from the observers side, it has a cavity in the 
SMBH direction. As a result, for every small BH in the 
volume ~ ttR 2 w R c \ there is amoment of fireball evolution 
when the rays, normal to the fireball surface, go through 
the Einstein radius of the BH. Therefore multiple lens- 
ing events are possible, like "mesolensings" in globular 
clusters [3]. But in contemporary models of GRBs the 
main flux of gamma radiation comes from jets, therefore 
the fireball surface is a subdominant source. In the case 
of A fireball scattered on the SMBH, different parts of 
its surface can be lensed similarly to [3l[- These lens- 
ings can explain the observed fast variability of the light 
curve of Swift J164449. 3+573451 during the first days. A 
characteristic feature of lensing is synchronous radiation 
in different bands. 



D. Quasi-steady radiation 

Multiple NS and BH collisions first produce dense flows 
of debris and then an extended hot gas cloud in the cen- 
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tral region of the cluster. Hot gas from the surrounding 
cloud falls onto the SMBH and radiates with the Ed- 
dington luminosity L+ — 47rGMsMBHTipc/cTT, where m p 
is the proton mass and <tt is the Thompson cross-section. 
This luminosity is 1.5 x 10 45 erg s _1 in the case of 10 7 M Q 
SMBH which is the minimum inter-peak luminosity dur- 
ing the first several days of Swift J164449. 3+573451. Ra- 
diation from slow gas accretion onto the newly formed 
SMBH requires the gas envelope to be thin for radiation. 



E. Multiple fireballs in the dense cloud 

The cluster of NSs and BHs is enveloped by a gas 
cloud originating from destruction of ordinary stars and 
partly from primordial gas in the galactic center. Mul- 
tiple GRBs are generated in the medium "prepared" by 
the preceding GRBs, so the evolution of multiple fireballs 
must be considered as in (32|, [33[. 

Collisions of compact objects produce the relativis- 
tic shock with Lorentz factor r 3> 1, which propagates 
ahead of the fireball pushing the gas in the nucleus. 
The shock becomes non-relativistic at the Sedov length 
l s = [3i?o/(47rp e c 2 )] 1 / 3 , where p e is the density of the gas 
in the envelope. The length l s gives an estimate for cav- 
ity radius, produced by the single fireball in the absence 
of the central SMBH. The multiple expanding fireballs in 
the cavity have the shape of thin shells and are separated 
by the distance R c — ci cap , where i cap corresponds only 
to NS-NS and NS-BH collisions in this case. The gas be- 
tween each two fireballs is swept up and compressed by 
the preceding one. Multiple fireballs produce the station- 
ary cavity |32| with radius i? cav = [N c Eq/ (Airpf-c 3 )] 1 / 2 , 
where c s is the sound speed in the gas. For the relativistic 
fireballs the gravitational attraction of the central SMBH 
is not important except for the later non-relativistic stage 
with r < 1. The SMBH causes a deformation of the 
fireball, opening the conditions for gravitational lensing. 
One can assume an equipartition of the magnetic field in- 
duced by turbulence and by the dynamo mechanism. In 
this situation the Fermi II acceleration mechanism oper- 
ates, and high energy particle creation is possible. The 
signatures of these processes are expected from the di- 
rection of Swift J164449.3+573451. 



V. CONCLUSION 

We propose a scenario of multiband powerful signal 
generation in an evolved star cluster in a galactic nu- 
cleus during the gravitational collapse of its dense cen- 
tral core. This model naturally provides necessary condi- 
tions for the event Swift J164449. 3+573451 because the 
huge power, duration, and variabilitymay indicate a more 
catastrophic event (or sequence of events) than tidal de- 
struction of a single star. 

In the typical case, the dynamical evolution of a central 
star cluster results in collisions and destruction of a sub- 
stantial fraction of stars and formation of a gas envelope 
[l2| . Simultaneously, a compact cluster of NSs and BHs 
can form in its evolution, due to sinking of massive stars 
to the center of the cluster and their explosions as super- 
novae [9l-[Il1] . The evolution of such a cluster leads to an 
avalanche contraction and gravitational collapse of the 
central region [l9j . This mechanism cannot explain the 
appearance of all SMBHs, especially in quasars at high 
redshifts, because of the lack of dynamical time. But 
star cluster collapses are nevertheless inevitable in many 
galaxies at low redshifts. Recently, the sub-dominant role 
of galaxies mergers in triggering the AGN activity was re- 
vealed [34J . Avalanche contractions and collapses of star 
clusters could be an alternative mechanism for AGN ig- 
nition at small redshifts along with tidal disruptions of 
stars or disk instabilities. 

Multiple dense gas releases and fireballs in NS-BH col- 
lisions can produce powerful peaks in the light curve of 
Swift J164449.3+573451. The additional inter-peak radi- 
ation comes from gas accretion onto the forming SMBH. 
We can expect that the properties of fireballs and gamma 
radiation in this model are distorted by their scattering 
on the growing central SMBH (the characteristic time 
scale is 2GM/c 3 ~ 100 s), and by gravitational lensing 
on compact objects in the cluster. 
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